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Abstract. We present a unified explanation for the observed properties of the isolated neutron star RX J0720.4-3125 by ob- 
taining a self-consistent model that accounts simultaneously for the observed X-ray spectrum and optical excess, the pulsed 
fraction, the observed spectral feature around 0.3 keV, and the long-term spectral evolution. We show that all observed proper- 
ties are consistent with a normal neutron star with a proper radius of about 12 km, a temperature at the magnetic pole of about 
100 eV and a magnetic field strength of 2 x 10'^^ G, value inferred from the observed period decay. The high magnetic field pro- 
duces a strong anisotropy in the surface temperature distribution. The observed variability of the effective temperature, strength 
of the spectral feature, and pulsed fraction are in good agreement with the predictions of our model in which the star is subject 
to free precession, producing changes in the angle between the magnetic field and the rotation axis of tens of degrees with a 
periodicity of 7 years, as pointed out by other authors on the basis of phenomenological models. In addition to the evidence 
of internal toroidal components, we also find strong evidence of non-dipolar magnetic fields, since all spectral properties are 
better reproduced with models with strong quadrupolar components. 
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1. Introduction 

RX J0720.4-3125 belongs to the family of radio-quiet iso- 
lated neutron stars (INS), a puzzling group of compact objects 
^ that during the last decade have provoked speculations about 
■ their real nature (neutron stars or strange stars) and forced us 
' to reconsider the thermal emission mecha nisms. RX J0720.4- 
3125 was discovered with ROSAT ( Haberletal. 1997>) . and 
its X-ray spectrum was soon found to be well described by 
a Planckian spectrum with temperature kT ~ 82 eV (k be- 
ing the Boltzmann constant). Similarly to the rest of INS, 
it is a nearby object (» 300 pc, Kaplan et al. 2003) and 
shows low interstellar absorption («// = 1-1.5 x 10^" cm"^). 
More interesting ly, it is a confirme d X-ray pulsar with a pe- 
riod of 8.391 s dHaberletal. 19971) and it is one of the two 
INSs with a r eliable measure of the peri od derivative P = 
7 X 10 ''^ss ' dKapIan & van Kerkwiik 2005 a,) . which implies 
a magnetic field of about B = 2.4 x lO'^ G Another impor- 
tant observational property, common to other INS, is that the 
observed optical flux is larger (about a factor of 6) than the 
extrapolation to the optical band of the best blackbody (BB) 
fit to the X-ray emission. This apparent optical excess flux of 



' The other object with a measure of the period derivative is 
RBS 1223, with P= 1.12x 10'" ^.?'', which implies B = 3.4x10" G 
teaplan & van Kerkwiik 2005bh 



INSs, first observed in RX J1865-3754 can be explained with 
the existence of large temperature anisotropies over the surface 
kons et al. 2002). In the case of RX J0720.4-3125, the evi- 
dence of anisotropic temperature distributions is strongly sup- 
ported by the observed X-ray pulsations with relatively large 
amplitude (~ 11%). 

More recently, the story of INSs has suffered a new 
twist when observations with XMM-Newton have revealed 
deviations from a pure BB spectrum in the form of absorption 
features observed in the 0.1 - 1.0 keV band. In the particular 
case of RX J0720.4-3125, a phase depen dent absorption lin e 
around 270 eV has been recently reported ("Habe rl et al. 2004 . 
This feature is normally associated with proton cyclotron 
resonant absor ption and/or bound-bound transitions in H 
or H-like He ('Haberl et al. 2003*, Ivan Kerkwiik et al. 2004h . 
Both require a magnetic field of ~ 5 x 10^^ G, consistent 
within with the dipole breaking estimate of 2.4 x 10'^ G, 
but only within a factor 2. Moreover, this estimate only 
reflects the large scale magnetic field structure, whereas the 
presence of smaller scale inner magnetic fields (i.e. strong 
toroidal components) seems very likely because MHD core 
collapse simulations show that toroidal magnetic fields are 
quickly generated by differential rotation and early stage 
convective motions can also create toroidal fields. Strong 
magnetic fields can afi'ect the emission properties of the 
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neutron star surface in multiple ways. For example, they can 
induce a phase transition tur ning the gaseous atmosphere 
into a liquid or solid state jLai 200 ih . which results in a 
spectrum that, besides a reduction in the emissivity compared 
to a BB, shows some particular features that could explain 
the observations. This scenario has been extensively studied 
by different authors jBrinkman n 1980>. 'Turolla et al. 2004 
IPerez-Azorm et al. 200a Ivan Adelsberg et al. 2005). In 
addition, if the magnetic field is high enough to induce 
the condensation of the atmosphere, it will also lead to 
very large anisotropics on th e surface temperature dis- 
tribution jGeppert et al. 2004t IPerez-Azorm et al. 20061 
iGeppert et al. 20061) . providing an attractive scenario to nat- 
urally explain the observed large optical excess and pulsed 
fraction of some isolated neutron stars such as RX J0720.4- 
3125 and RBS 1223. The case of RX J0720.4-3125 is also 
particularly interestin g because of its clea r spectral evolution 
jde Vries et al. 2004t iHaberl et al. 2004t IVink et al. 20041) 
which has been recently associated with a ~7 yr precession 
period of the neutron star (Haberl et al. 2006). 

The paper is organized as follows. In Section 2 we review 
the main properties of our theoretical model (magnetic field 
configuration, temperature anisotropy, etc.). In Section 3 we 
summarize the X-ray observations used in this work and in 
Section 4 we revisit the phenomenological models previously 
adopted to describe the data. In section 5, we discuss and mo- 
tivate the choices made to limit the parameter-space of our 
models, and we present our main results, obtained by apply- 
ing our realistic and physically motivated models to the avail- 
able X-ray data. In Section 6 and 7 we discuss the evidence of 
precession and the interpretation of the excess optical flux, re- 
spectively. Finally, in Section 8, we summarize our results and 
discuss open questions and caveats. 

2. The theoretical model. 

In previous papers jPerez-Azorm et al. 2005[ 

IPerez-Azorm et al. 20061) we have presented the results 
of detailed calculations of the temperature distribution in the 
crust and condensed envelopes of neutron stars in the presence 
of strong magnetic fields, by obtaining axisymmetric, station- 
ary solutions of the heat diffusion equation with anisotropic 
thermal conductivities. Having explored a variety of magnetic 
field strengths and configurations, we concluded that variations 
in the surface temperature of factors 2-10 are easily obtained 
with B X lO'^-lO'"* G whereas the average luminosity (and 
therefore the inferred effective temperature) depends only 
weakly on the strength of the magnetic field. 

Nevertheless, the luminosity is drastically affected by the 
field geometry, in particular by the existence of a toroidal com- 
ponent. Moreover, if the magnetic field is strong enough to in- 
duce the condensation of the surface, the condensed surface 
models also predict the existence of a spectral edge that for 
B ^ lO^^-lO''* falls in the range 0.2-0.5 keV and that can be 
consistent with an absorption line such as that reported in RX 
J 0720.4-3125 ( iHaberl et al. 2004. We refer to the in terested 
reader to the previous work jPerez-Azorm et al. 2006*) for de- 
tails about the calculations, and we sketch here only the main 



features of the magnetic field geometry and surface tempera- 
ture distribution. 

The structure of relativistic stars with both poloidal and 
toroidal magnetic field components has recently been studied 
(loka & Sasaki 2004). It was shown that all quantities can be 
determined from a stream function that satisfies the relativis- 
tic Grad-Shafranov equation. In the linear regime (weak mag- 
netic field, in the sense that deformations are small), the Grad- 
Shafranov equation becomes simpler and, under the assump- 
tions of axisymmetry and stationarity, the general interior solu- 
tion of the magnetic field has the form (loka & Sasaki 2004h 

/ cose smdd'¥{r) sind \ 

B = Bo2^^(r), ^'{r)\ (1) 

\ r or r ) 

where yu is a constant with the interpretation of the wavenumber 
of the magnetic field. In general, the /-component of the stream 
function must satisfy the following differential equation 

^^2 +1'" I ^('■^ = P^O ' 

where p is the energy density, ao = ao(r) is a function that 
depends on the boundary conditions and we have omitted the 
relativistic factors for clarity. In the case ju = the purely 
poloi dal cas e is recovered, as discussed for example in Konno 
et al. d 19991) for the weak magnetic fi eld lim it, and more thor- 
oughly in Bonazzola & Gourgoulhon (1996) for the fully non- 
linear case. For the general case with both poloidal and toroidal 
components, a similar non-linear (although New tonian) analy- 
sis is presented in Tomimura & Eriguchi (120051) where Eq. ^ 
is generalized by considering fi - iu{x) as a general function 
of jc = r sin 0A^, where is the component of the vector 
potential. 

It is not the purpose of this paper to review all previous 
calculations about the structure of neutron stars with magnetic 
fields in detail, but we must stress the more relevant issues that 
lead us to chose a particular magnetic field structure. The key 
point is that under the assumptions of stationarity and axisym- 
metry the electric cuiTent (source of the magnetic field) can 
only have two components, namely a purely toroidal current 
which produces a poloidal field, and a second term propor- 
tional to the magnetic field, and therefore resulting in a force- 
free field with both poloidal and toroidal components. Thus, the 
magnetic field is a combination of a force-free component and 
a purely poloidal one. Since the deformation of neutron stars 
with the magnetic fields intensity of interest (lO'^-lO'"^ G) are 
known to be fairly small, it is reasonable to assume that the 
field structure can be approximated by the force-free solution 
(that is flo = in Eq. ([3). In this case (and in the Newtonian 
limit) Eq. (|2} is a form of the Riccati-Bessel eqti ation, which 
has analytical solutions ( Perez-Azorin et al. 20061) . 

In this paper, we will consider force-free solutions with 
fi chosen to confine the magnetic field to the crust and outer 
regio ns, as some of the models discussed in previous works 
(Geppert et al. 2004 I Perez-Azorfn et al. 20061) . where it was 
shown that this crustal confined configurations lead to large 
anisotropics in the surface temperature distribution. In particu- 
lar, for our neutron star models, we have considered two cases 
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f(0) = /3rfCose+/?^(3cos2 6i- l)/2^ 



(5) 
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Fig. 1. Surface temperature profiles as a function of the polar 
angle for different magnetic field configurations according to 
Eq. |5l We have taken = 1.34 km"' and fj.2 - 3.877 km"'. 
The different lines correspond to different relative strengths 
between the purely dipolar (J3^ = 1) and purely quadrupolar 
(/3d = 0) models. 

in which yU - 1.34 km"' and fj. = 3.87 km"' depending on the 
considered magnetic field configuration, either purely dipolar 
or quadrupole-dominated (see Perez-Azorm et al. (2006) for 
plots of the magnetic field lines). Notice that jiR is roughly the 
ratio of the toroidal to the poloidal field, while, as mentioned 
above, //"' has the interpretation of the typical length scale of 
the magnetic field. 

It turns out that for models in which the toroidal magnetic 
field covers the whole crustal region (as in force-free configu- 
rations) most of the temperature gradient is located in the en- 
velope. Therefore, the classical semi-anal ytic te mperature dis 
tribution derived by Greenstein & Hartke 



ytic te n 
il983h 



T^, Jcos^ ^« + ~ ^^"^ 



(3) 



remains a good approximation. Here Ob is the angle between 
the magnetic field and the normal to the surface. A simple but 
accurate analytical approximation to the temperature distribu- 
tion for dipolar force-free models, which is more realistic that 
discontinuous two temperature models, is 



A A COS A 

= 2a ^..1 ■ + 

cos-^ + ^ sm- 9 



(4) 



is the 



where Tp is the temperature at the magnetic pole, T„ 
minimum temperature reached in the surface of the star (typ- 
ically Tmin < lOFp) and the parameter a takes into account 
the relative strength between poloidal and toroidal components. 
For a dipolar field a = 0, while for force-free models with mag- 
netic field confined to the crust and envelope a ^ 250. An anal- 
ogous formula when quadrupolar components are added can 
also be derived: 

.4 



T\0) 



'f(0y + 2±i(j3d+/3gcos6y- sin^e 



+ TZ 



where fid and fiq - \ - fid are the relative weights between the 
dipolar and quadrupolar components. We have found that this 
approximate expressions are accurate to less than 3%. 

It must be mentioned that this force-free structures 
can be continuously extended into a twiste d magnetosphere 
jLvutikov et al. 2002t iLvutikov & Gavriil 20 06). which has 
been proposed to be the explanation of the infrared and optical 
emission of magnetars (Durant & van Kerkwiik 2005). Within 
this model, the twisted magnetic field is present not only in the 
magnetosphere, but also in the envelope and probably in the 
crust of neutron stars. 

As an example of temperature distributions, size of the 
hot spots, and effect of higher order multipoles, in Fig. ^ 
we show surface temperature profiles as a function of the 
polar angle for different magnetic field configurations vary- 
ing from the purely dipolar (fid - 1) to purely quadrupolar 
(fid - 0) models. The main features are: i) purely dipolar 
models show nearly equal antipodal hot polar caps with only 
minor differences due to the Hall term in the diffusion equa- 
tion (Geppert et al. 2006: Perez-Az orm et al. 2006), ii) when a 
quadrupolar component is present, the north/south asymmetry 
is broken by the different temperatures and sizes of the polar 
caps, and iii) a dominant quadrupolar component implies also 
the existence of a hot equatorial belt, that can be displaced from 
the equator depending on the relative strength of the dipolar 
and quadrupolar components. Of course similar but more com- 
plex geometries could be obtained by including higher order 
multipoles. In general, smaller values of ji correspond to larger 
angular sizes of polar caps and/or equatorial belts. 

3. Summary of XMM-Newton observations 

RX J 0720.4-3125 has been observed several times by XMM- 
Newton and we focus here on EPIC data (pn and MOS 1) from 
the 8 publicly available observations. The XMM-Newton ob- 
servations span a period of about 5.5 years starting from May 
2000 (Rev. 78) to November 2005 (Rev. 1086). Here we ig- 
nore the observations performed with the cameras operated in 
Small Window because of the worse calibration with respect to 
Full Frame science modes. The pn data of Rev. 078 are also ex- 
cluded from the analysis due to problems in the SAS 6.5.0 re- 
duction pipeline (see http : //xmm. vilspa . esa.es/ for de- 
tails). The remaining data provide a homogeneous set and have 
all been collected with the cameras operated in Full Frame 
mode with the Thin or Medium filter applied and a summary 
of the observations is presented in Table ^ The MOS 2 data 
are consistent with MOS 1 and do not add relevant information 
to our analysis. Photon pile-up has been minimised by follow- 
ing standard procedures (see e.g. Haberl et al. 2004). Since the 
softest energies suffer from calibration uncertainties, we con- 
sider the 0.18-1 .2 keV band only, after having checked that the 
inclusion of data down to 0. 13 keV does not change our results 
in any noticeable way (worsening the statistics in a similar way 
for any spectral model). 

In the following we present X-ray data-analysis for all 
observations by using three different sets of models, a phe- 
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nomenological one (following previous works by e.g. Haberl 
et al. 2004; 2006) and two realistic models obtained through 
detailed numerical simulations in the theoretical framework de- 
scribed above. In all three cases we also consider a sub-set of 
the pn-FF observations in which we force the absorbing col- 
umn density to be the same in all observations. For the two re- 
alistic models, we also force the model normalization (related 
to the NS size and distance) and the magnetic field intensity to 
be the same in all observations. 

4. Phenomenological fits revised 

As already shown in previous w orks dde Vries et al. 2004t 
iHaberl et al. 2004t IVink et al. 2004 . the X-ray spectral shape 
of RX J 0720.4-3125 is not that of a pure absorbed black- 
body (BB). Deviations from a thermal spectrum are seen in 
the 0.2-0.6 keV band and they are well fitted by models in 
which a Gaussian absorption line is added to the pure (ab- 
sorbed) thermal spectrum. The line was interpreted as cy- 
clotron resonance s cattering of protons in the NS magnetic field 
llHaberl et al. 2004 ). While the above interpretation is certainly 
allowed by the data, the exact shape of the absorption struc- 
ture cannot be constrained with high confidence from the X-ray 

J ' 

data (IVink et al. 20 04) and possibly describes in a phenomeno- 
logical way the deviations from a pure thermal spectrum. We 
point out here that BB plus absorption line fits to the X-ray 
data are phenomenological almost by definition since they can- 
not account simultaneously for the observed Optical and UV 
data, underestimating the Optical flux by a large factor. This is 
the reason why in the following we shall call this description 
of the X-ray data "phenomenological". 

We repeated this analysis by considering all phase- 
averaged XMM-Newton observations and by describing the 
spectra with a simple BB plus Gaussian absorption line model, 
photo-electrically absorbed by a column density of gas, and 
our results are reported in Table |2l Errors are given at the 90% 
level for one interesting parameter The Gaussian absorption 
line width cannot be constrained with high confidence and it 
is fixed to its best-fit average value to all observations (75 eV, 
similar to the 64 eV width imposed by Haberl et al. 2004). We 
point out that if the Gaussian width is let free to vary, the ab- 
sorption line energy becomes unconstrained in many cases. All 
other parameters are free to vary. The final statistics is good 
for both the pn and MOS 1 cameras and the absorption line is 
required at more than the 99.99 % confidence level in all ob- 
servations except Rev. 078, 175, and 622. 

As mentioned, we also consider a smaller but more homo- 
geneous sub-set of data from the EPIC-pn camera only, always 
operated in Full Frame. In this case, since inter-calibration 
between different instruments is not anymore a concern, we 
force the absorbing column density to be the same in all ob- 
servations, as it seems reasonable. Our results are presented in 
Table Island suggest that both the BB temperature and the ab- 
sorption line equivalent width increase with time. Both quanti- 
ties increase up to Rev. 815 and decrease in the last observ ation 
(Rev. 1086), confirming previous results ( iHaberl et al. 20 06). 
The temperature variation (from ~ 85 eV to ~ 91 eV) may 
indicate that the hot spot(s) responsible for the X-ray pulsa- 



tions of the source has a different viewing angle (possibly due 
to precession), yielding to higher and hig her effective temper- 
atures (as proposed e.g. in lVink et al. 20041) . Evidence for pre- 
cession in R X J0720.4-3125 ha s been indeed strengthen in a 
recent work dHaherl etal. 20061) where a ~7 yr periodic varia- 
tion of the BB temperature has been claimed, based on similar 
fits (BB and Gaussian absorption line) to the X-ray data. Our 
results on the variation of the BB temperature are consistent 
with the inferred period. 

5. Spectral analysis with realistic, self-consistent 
models. 

In the previous Section, we discussed simple BB plus absorp- 
tion line fits to the X-ray spectra of RX J0720.4-3125. We 
remind here that though such models do provide a very good 
description of the X-ray data, they cannot explain in a self- 
consistent way the Optical flux which is observed to exceed by 
about a factor 6 the predicted one. In the following we present 
an attempt to describe the X-ray data with synthetic spectra ob- 
tained through detailed numerical simulations which can also 
simultaneously provide a good description of the Optical/UV 
data. 

However, the task of fitting observational data in a multi- 
parameter space is numerically time-consuming (because it re- 
quires computing a large number of models to cover a sufficient 
number of grid points) and theoretically complex (because it is 
difficult to avoid falling in local minima in a highly dimen- 
sional parameter space). For this reason we have not attempted 
to solve the problem using brute force but rather to discrimi- 
nate which parameters are more relevant for each observational 
fact. We have considered a fiducial NS model fixing the mass at 
1 .4 Mq and the radius at 12.27 km. While the first choice is well 
motivated by observations of NSs in binary systems the radius 
of NSs is much less certain as revealed by our ignorance of the 
details of the inner nuclear structure and equation of state. 

In the present case, however, we are interested in the X-ray 
(and Optical) emission from the star surface which, in absence 
of spectral features that can reveal the redshift at the surface, 
will depend very little on the assumed radius. Indeed, our mod- 
els are not significantly affected if the NS radius is varied in a 
sensible range (e.g. 10-15 km). This motivates and justifies the 
choice of a fiducial model which greatly simplifies the numeri- 
cal task of computing large grids of models if all parameters are 
allowed to vary. The magnetic field intensity is defined at the 
pole and we consider variations in the range of 0.5 - 6 x 10'^ G, 
as suggested by the value of 2.4 x lO'^ G estimated from the 
period decay. As for the magnetic field configuration, in the fol- 
lowing we wiU explore two different cases, both confined in the 
outermost regions (crust and envelope): a purely dipolar mag- 
netic field with ji - 1.34 km and a quadrupole-dominated 
one - 0.05) with yu - 3.87 km"'. This is only a small sub- 
set of the possible configurations which has been chosen after a 
relatively time-consuming assessment of which configuration 
provides a better description of the X-ray data. In particular, 
(very) different values for the radial length scale of the mag- 
netic field (ju) resulted in inconsistent results in which the ra- 
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Fig. 2. Contour plots for a purely dipolar magnetic field, with 
HH = 1.33 X 10^0 cm-2, = 115 eV and Bp = 1.8 x 10'^ G, 
of the pulsed fraction for a realistic surface temperature distri- 
bution. S is the angle between rotation and magnetic axis, and 
O the angle between the rotation axis and the observer direc- 
tion. The flux has been obtained by integration over the whole 
surface taking into account General Relativistic light bending 
effects. 

dius of the NS, as inferred from the X-ray spectral fits, was 
different from the fiducial model one. 

In the future, we plan to build a much more extensive set 
of grid models which will be available for X-ray spectral fit- 
ting purposes with most of the parameters free to vary. This 
will result inevitably in a better statistical description of the X- 
ray data than possible at present. We consider the present work 
as a first attempt to reproduce the available observational as- 
pects with a realistic self-consistent model, without claiming 
that the proposed model(s) is (are) the best-possible solution. 
However, as discussed later, the results inferred from our anal- 
ysis provide several insights and suggest the direction to take 
with the goal of improving further the theoretical models and 
our understanding of INSs in general. 

Having set up a baseline model, we have built tabular 
XSPEC models ^ as a function of the pole temperature and 
magnetic field strength for several orientations within the range 
allowed by the pulsation profiles, as discussed below. We then 
consider all the public available XMM-Newton observations 
of RX J 0720.4-3125 and compare our realistic (and limited 
in parameter-space) models to the data and to the more phe- 
nomenological description of the X-ray spectra discussed in 
the previous Section. 

5.1. Purely dipolar magnetic fields. 

We begin our analysis with a purely dipolar magnetic field con- 
figuration. This is mainly suggested by the observed regular si- 
nusoidal shape of the light curve of RX J0720.4-3125 (when 
folded at the spin period of 8.391 s). In the framework of our 
models, the spectral shape and time-evolution mainly depends 

^ Tabular models are available upon request to the authors 



on the relative orientation between the magnetic axis, the ro- 
tation axis, and the observer line of sight. These can be con- 
strained by using the observed pulsation profile and amplitude 
of RX J0720.4-3125 as detailed below. In Fig.|2lwe show con- 
tour plots of the pulsed fraction, for one of our calculations, as 
a function of the angle between rotation and magnetic axis (S) 
and between the rotation and observer direction (O). 

First, we can classify the light curves in two groups. The 
models with large O + S are characterized by a non sinusoidal 
pulsation, or even two visible maxima in the pulse profile; this 
is because we actually can see both poles in each period if 
O + S > 90°. Therefore, INSs which exhibit non-sinusoidal 
pulse profiles do lie in this region. On the other hand, the mod- 
els with small O + S always show one single peak in the pulse 
profile which is very close to sinusoidal when either one of the 
angles is small. This is probably the case of RX J0720.4-3125, 
that shows a very regular pulsation profile. 

Next, we can reduce the range of angles to those consistent 
with the observed pulsed fraction. For a nearly spherical neu- 
tron star, the rotation axis is essentially aligned with the (vec- 
tor) total angular momentum of the star, which is conserved. 
Therefore, the variation of the angle O with time is too small 
to be observable, but the star wobbles around its symmetry 
axis (in general different from both, rotation and magnetic axis) 
with a free precessio n timescale of a few vear s for oblateness 
of the order of 10"^ (I.Tones & Andersson 200 ll) . We have tried 
different orientations in the range of O and S allowed by the 
observed pulsed fraction. We find that the observed pulsed frac- 
tion values are best reproduced if (9 = 12° and S varies in the 
range 30°-60° (i.e. 180° -S in the range 120°-150°). With this 
choice of parameters, the pulsed fraction of our models always 
lies in the range between 9% and 13%, consistent with the ob- 
served values. Therefore, we will focus on the vertical shaded 
region, fixing O and allowing for variations in S. We have build 
tabular XSPEC spectral models by fixing O - 12° with the 
following free parameters: pole temperature KT^ok, magnetic 
field intensity Bp, magnetic field orientation S, and model nor- 
malization 7?oo/D3oo. The model is applied to all available data- 
sets with the addition of photo-electric absorption. 

As for the phenomenological model described above, we 
first use all the pn and MOS 1 observations with all parameters 
free to vary. Our results are presented in Table |3 The (unred- 
shifted) polar temperature varies approximately from 110 eV 
to 120 eV about 25% larger than the BB temperature of Table 
|2j while the average effective temperature over the NS surface 
is about 35-41 eV, more than a factor 2 smaller than the tem- 
perature inferred from BB fits. This is due to the very large 
anisotropy over the surface temperature (a factor 10 between 
pole and equator) induced by the magnetic field. The X-ray 
spectrum is dominated by the small hot polar area, while the 
extended cooler equatorial belt is responsible for most of the 
optical flux which is well reproduced (see Section 5.5). 

The statistical quality of the fits obtained with the realistic 
models is comparable to that obtained with BB plus Gaussian 
models (see Table |2] for comparison). However, the accuracy 
on the most relevant parameters (S, KTpoif^, and Bp) is limited 
by the fact that the inferred variations of ne drive the fitting 
results, strongly affecting the soft energy band where most of 
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the photons are collected. Given the high proper motion of the 
source a variation in the absorbing column cannot be excluded, 
but the observed variation seem to occur almost randomly and 
the pn and MOS 1 data often give inconsistent results for the 
same observation. It is much more realistic to assume that the 
column density is the same in all observations. 

To overcome the inter-calibration uncertainties, we con- 
sider the pn-FF observations only. As for the BB plus Gaussian 
fits described above, we force the absorbing column density to 
be the same in all observations. We also force the model nor- 
malization (directly related to the NS radius and distance) and 
the magnetic field intensity to be the same. Our results are pre- 
sented in Table |5] As it can be seen by comparing with the 
results in Tabled the statistical quality of the fits is worse than 
in the BB plus Gaussian case. The best-fitting magnetic field 
intensity Bp is of the order of 1.4 x lO'^ G, about a factor 2 
smaller than that inferred from the observed period decay, pos- 
sibly indicating that our model is not an extremely accurate 
description of the magnetic field geometry. 

Despite the above problem, a clear long-term evolution of 
the magnetic field orientation S emerges from our results. S in- 
creases with time up to Rev. 815, while it decreases in the last 
observation (Rev. 1086). By considering the vertical shaded 
area in Fig.|2](or its equivalent complementary at kS), this so- 
lution suggests that the pulsed fraction must exhibit a long term 
evolution, with a maximum pulsed fraction around S ~ 55°, 
corresponding to the latest observations. We point out that the 
S long-term evolution has the same behaviour shown by the 
BB temperature and the absorption line EW (see Table |3}, and 
(see below) pulsed fraction. Our realistic model is thus trying 
to reproduce the observed variations with a change in the mag- 
netic field orientation, which is what is expected if precession 
is responsible for the long-term variability. 

The model explored so far seems to be able to describe the 
bulk of the observed long-term variability in terms of preces- 
sion and is therefore very promising. However, the quality of 
the fits is not very satisfactory and we thus explore in the fol- 
lowing another case in which a strong quadrupolar component 
is added to the magnetic field configuration. 

5.2. Quadrupole-dominated magnetic fields. 

In Fig. 13 we show contour plots of the pulsed fraction for 
one of the models with a dominant quadrupolar component 
(J3d = 0.05). The very different surface temperature distribution 
(two hot spots plus a hot belt close to the equator, as discussed 
in Section 2) produces significant differences with respect to 
the purely dipolar case. First, the allowed maximum pulsed 
fraction is much higher, up to 45%. Second, the north/south 
symmetry has been broken, and the results are different when 
changing either S or O hy n - S or n -O (the symmetry with re- 
spect the simultaneous interchange of both angles by their sup- 
plementary is kept). Third, a portion of the hot belt will now 
be visible, which will have consequences on the spectrum. We 
proceed as in the previous case, and we first localize the region 
in which we are allowed to vary the angles O and S, keeping the 
pulsed fraction at about 11%, and close to a sinusoidal shape. 




Fig. 3. Contour plots of the pulsed fraction for a realistic sur- 
face temperature distribution produced by a quadrupole domi- 
nated magnetic field (J3d = 0.05). The parameters of the model 
are Tp = 115 eV ,Bp = 1.8 x 10'^ G and = 1.2 x 10^" 
cm"^. The flux has been obtained by integration over the whole 
surface taking into account General Relativistic light bending 
effects. 



as observed. This corresponds to O - 11° and S = 20° - 65° 
(see Fig.E). As mentioned, O = 169° and S = 115° - 160° 
produces the same results. 

We proceed as above, first applying the model to all obser- 
vations and letting all parameters free to vary, and then consid- 
ering the sub-set of pn-FF data only, forcing the absorbing col- 
umn density, the model normalization, and the magnetic field 
intensity to be one and the same in all observations. Our re- 
sults are presented in Table |5] and Table ^respectively. When 
all parameters are free to vary, the statistical quality of the fit 
is comparable to that obtained with the phenomenological and 
purely dipolar-field models (Tables |2 and 0} but, as with the 
purely dipolar model, the variation of the absorbing column 
density is a possible sign of a degeneracy in the parameters. 

We thus regard the pn-FF observations in which the pa- 
rameters that should remain constant have been forced to be 
the same in all observations as the most reliable sub-set to in- 
fer information on the remaining model parameters (Table 0. 
In this case, we obtain a very significant improvement with re- 
spect to the purely-dipolar model and we approach the statis- 
tical quality of the phenomenological fits. With the addition of 
the quadrupolar component, the best-fitting magnetic field in- 
tensity raises to 1.8 + 0.2 x lO'"* G, much closer to the value of 
~ 2.4x 10'^ G inferred from the period decay. As for the dipolar 
model discussed above, the magnetic field orientation S clearly 
shows a long-term evolution similar to the n - S evolution of 
the dipolar model. The issue of the long-term evolution of the 
magnetic axis orientation will be discussed later in more detail 
in the framework of the precession interpretation for the vari- 
ability of RX J0720.4-3125. 
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5.3. The hardness ratio anti-correlation 

The better statistical quality of the fits to the X-ray data and 
the more consistent value of the magnetic field intensity are, up 
to some extent, an indication of the presence of higher order 
multipolar components. However, as recently discussed thor- 
oughly in (Zane & Turolla 2006) there is a more conclusive ob- 
servational fact that allows us to distinguish between the purely 
dipolar and the multipolar case. A striking feature of the X-ray 
light curves of RX J0720.4-3125 is the clear anti-correlation 
of the hardness ratio with the pulse profiles in both the hard 
and the soft band. This feature cannot be explained with ax- 
isymmetric models that have also a north/south symmetry. A 
possible explanation that has been proposed is a model with 
two hot spots with one of them displaced by about 20° from 
the south pole (Haberl et al. 2006). 

Motivated by the well established hardness ratio anti- 
correlation, we explored whether our models can reproduce the 
observed behaviour In Fig.l^we show the folded light curves 
in the hard and soft bands, and the hardness ratio for one of 
the quadrupolar-dominated models used in fitting the X-ray 
data. The observed anti-correlation of the hardness ratio with 
the pulse profiles is qualitatively very well reproduced by our 
model. We point out that such good agreement is impossible to 
obtain with pure dipolar fields in which the North/South sym- 
metry is not broken. We have not attempted at this stage to find 
the model that best fits the light curves and the phase averaged 
spectra at the same time, which is a formidable numerical task, 
especially when dealing with realistic models that (each one) 
need computational effort to be produced. However, the good 
qualitative agreement is indicative that we could be seeing a 
neutron star with a hot spot and a hot belt, as predicted by 
quadrupole-dominated models, as an alternative to the case of 
a displaced hot spot. 

If the interpretation that all variability is due to free preces- 
sion of the neutron star is correct, one should be able to find a 
unique model (temperature, magnetic field configuration, etc.) 
that explains all observations by only varying the relative ori- 
entation and, in our opinion, this model approaches that solu- 
tion. Our results point in the direction that most of the variation 
is actually explained by precession of the neutron star. Notice 
that a precession timescale of a few years has also already been 
reported for some pulsars (.Link & Epstein 2001.) . 

5.4. Precession 

As mentioned, evidence for precession in RX J0720.4-3125 
has been recently reported by Haberl et al. (2006). The evi- 
dence was based on results from BB plus absorption line fits to 
the available X-ray data in which the BB temperature shows a 
long-term evolution best described by a sine wave with a pe- 
riod of 7.1 + 0.5 yr. We consider our results for the dipolar and 
quadrupole-dominated models (see Table |5] and and try to 
fit the time evolution of the angle S which is the relevant quan- 
tity as far as precession is concerned. The results are shown 
in Fig. |5l Notice that the right panel shows n - S for better 
comparison. This is just a matter of the arbitrary choice of co- 
ordinate origin of the polar angle. Both are consistent with the 
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Fig. 4. Pulse profiles in two energy bands (soft:0. 12-0.4 keV, 
hard:0.4-1.0 keV) for a reahstic model with O ^ 1 1, S = 54, 
Tj, = 115 eV ,Bp = 1.8 X 10^^ G, rin ^ 1.2 x lO^^ cm^^ and 
/3d = 0.05. The third panel is the hardness ratio (HARD/SOFT). 
The bottom panel shows the observational data corresponding 
to observation rev. 534. Notice that we have not attempted to 
perform a fit of the pulse profiles, it is simply a comparison 
with the model obtained from spectral fitting. In this figure we 
have taken into account the response matrix of the instrument 
and the absorption. 

interpretation of a precessing neutron star with a precession pe- 
riod of 7 years and a relatively large wobbling angle of ^ 20°. 
Both the dipolar and quadrupole-dominated cases are very well 
described by a sine wave but, as mentioned earlier, the hard- 
ness ratio anti-correlation strongly favors model with a signif- 
icant quadrupole component breaking the north/south symme- 
try. Our results are consistent with the conclusions drawn in 
llHaberl et al. 20061) with phenomenological models. 

5.5. Pulsed fraction long-term evolution 

If precession is the right interpretation fro the long-term spec- 
tral evolution of RX J0720.4-3125, the observed pulsed frac- 
tion should also show long timescale variability associated 
with the variation of the hot spot(s) effective area. We have 
then extracted 0.12-1.2 keV light curves for the pn observa- 
tions used above (the Full Frame ones), computed the period 
of the observed pulsations (always consistent with 8.391 s), 
and obtained the pulsed fraction for each observation. Our re- 
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ata (see text in figure). 



suits suggest that the pulsed fraction is not constant in time, 
but increased from 10% to 12% on the few years timescale 
spanned by the observations. We measure a pulsed fraction of 
10.1 + 0.5% in Rev. 175, 11.1 + 0.4% in Rev. 533 and 534, 
and 11.7 + 0.5% in Rev. 815 and 11.9 + 0.5% in Rev. 1086. 
The beahaviour is suggestive of a correlation between the BB 
temperature, the absorption line equivalent width (or deviations 
from the BB spectrum) and the pulsed fraction, as expected in 
the precession model. 

In Fig. |6l we show the pulsed fraction long-term evo- 
lution with t he addition of a first data point from Rev 78 
(iHaberl et al. 2004). The solid line represents a fit with a sine 
wave in which the period is fixed at 7 yr, while the amplitude 
is a free parameter. We measure a modulation amplitude of 
1 .0 ± 0.3%. Although the quality of the data is not high enough 
to claim a periodic evolution of the pulsed fraction, the obvi- 
ous consistency with a period of 7 yr (see Fig.|6j strongly sup- 
ports the precession model. We point out that the pulsed frac- 
tion long-term evolution is independent of the adopted spec- 
tral model and of the instrument calibrations and represents 
thereby an independent indication that precession is at work 
inRX J0720.4-3125. 



6. Optical flux 

An issue under debate in the condensed surface models is 
related to the optical flux. When the eff'ect of motion of 
the ions is neglected the optical flux is very ni uch de- 
pressed (Tur olla et al. 2 004: Per ez-Azorm et al . 2005), but ig- 
noring the motion of ions at low energies is not justified. In fact, 
a simple treatment of the motion of the ions as free particles 
leads to very different results (Ginzburg, 1970). This simple 
approach (to consider them as free particles) could be approx- 
imately correct if the surface is in a liquid state, but it is not 
proper if matter is in solid state and ions are placed in a lattice. 
However, considering them as totally fixed to the lattice can be 
as wrong as letting them to move freely. In reality, ions do ac- 
tually move and resonant plasmon or phonon excitations can 
be important. Moreover, the strong magnetic fields can modify 



the motion of the ions (the ion cyclotron radius is much smaller 
than the separation between ions in the lattice). 

The usual Drude model that describes the microscopic re- 
sponse to electromagnetic fields in metals predicts that the di- 
electric function e at low energy behaves as 



(6) 



(jj{ti) + iT) 



where ojp is the electron plasma frequency (w^ = Ane^nc/me) 
and r is a relaxation coefficients. From this formula it is clear 
why at very low energy (cj « cjp) the absorption is strongly 
suppressed. In the other limit, assuming the motion of the ions 
is free because they are not bound in a lattice, the effect of 
the magnetic field is to modify the dielectric function (very 
schematically, and neglecting damping) as follows 



1 



(7) 



where (. and ajg i are the electron and ion cyclotron frequen- 
cies. Therefore, when o) < ojbj the mode which had a large 
imaginary part of the refra ction index is not damped any more 
jPerez-Azorin et al. 20061) . 

But real materials are known to behave in a more com- 
plex way than this. The simplest extension of the model that 
can account for collective excitations of the lattice at charac- 
teristic frequencies (wqj) is called Drude-Lorentz model, and 
it consists in the combination of a Drude term with a number 
of Lorentz oscillators 



e= 1 - 



+ 2- 



(8) 



It is clearly out of the scope of this paper to discuss the 
microphysical behaviour of this condensed (liquid or solid) 
atmospheres, and we do not know much about the charac- 
teristic frequencies of the lattice and how this couples to the 
magnetic field. We just want to point out that several impor- 
tant energy scales; the ion cyclotron energy hu)B,i, kT, and 
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Fig. 7. Spectral energy distribution of the best fits of the revo- 
lution 534 for phenomenological model (Table|3j black solid), 
realistic quadrupolar dominated (Table0 blue solid), and real- 
istic dipole (Table|5j red solid). The phase-averaged spectra of 
both models are very similar, and it is necessary to analyze the 
spectral phase evolution to discriminate. The dotted line is the 
optical tail of the BB model. Solid and dashed lines in the op- 
tical band correspond to models with free and frozen ions, re- 
spectively. Available optical and UV data dKaplan et al. 20031) 
are also shown (crosses) for comparison. 



fi(^p,ion are similar (about 0.1 keV) under our typical con- 
ditions and, consequently, a careful analysis of the emissiv- 
ity in the optical band is needed before we can establish se- 
vere constraints on the models. As a compromise, since the 
free ion treatment predicts an emissivity in the optical band 
(a - I corresponds to a blackbody) of about a = 0.7 
and a frozen ions model predicts approximately a - 0.2 
jPer ez-Azor m et al. 20051 Fvan Adelsberg et al. 2005 ^. we will 
show both limits in the following figures, so that reality must 
be places in between this two different possibilities. 

In Fig. we show the unfolded spectra for the best fits 
to the X-ray spectrum of Rev. 534 together with the available 



optical-UV data d Kaplan et al. 20031) . We have not attempted 
to fit simultaneously the optical and X-ray data, but the good 
agreement is evident. For comparison, we also show (dotted 
lines) the best BB plus Gaussian fit for the Rev. 534 data to il- 
lustrate its similarity with the realistic model in the X-ray band 
and to highlight the well known problem of under-predicting 
the optical flux. 

Notice that the best fits to the optical observations at lower 
energies (Jul 2001) show a somewhat higher flux than the latter 
optical observations (Feb 2002). The fact the the optical obser- 
vations are not consistent with a pure Rayleig h-Jeans tail has 
been discussed in detail in dKaplan et al. 20031) who found that 
the best fit consists of two BBs plus a power law. However, 
they also comment that the spectrum can be consistent with 
a Rayleigh-Jeans tail if the deviations have a temporal nature, 
which is thought to be unlikely because the X-ray flux is con- 
stant. Interestingly, the fact that all observations fall in between 
the two lines predicted by the/ree ions frozen ions models 
may be telling us that more effort must be placed in understand- 
ing the physics of this emission processes. 

7. Summary and conclusions 

We have shown that realistic models of neutron stars with 
strong magnetic fields are consistent with the observed X-ray 
and optical spectra, the observed deviation from a pure ther- 
mal spectrum in the X-ray band, and the long term variability. 
Although we do not exclude the presence of a resonant proton 
cyclotron absorption or bound-bound transitions in H, we do 
not need to invoke it to explain the observed spectral edge. We 
stress that our model calculates the surface thermal emission, 
which could further modified by the interaction of the radiation 
with the plasma in the magnetosphere. 

Once the magnetic field configuration is given, we are 
able to obtain a self-consistently calculated thermal spectrum 
that reproduces reasonably well all available observations. Our 
analysis indicates that most of the long-term spectral varia- 
tion of the source can be explained in terms of neutron star 
precession. If our interpretation is correct, we predict that 
the pulsed fraction, effective temperature, and absorption line 
equivalent width as obtained with phenomenological models, 
will decrease in subsequent observations of the source, until 
reaching similar values to Rev. 175. This begins to be appar- 
ent in the recent observation Rev. 1086, made public while 
we were revising the final version of this manuscript. Similar 
conclusions have b een reached independently by other groups 
dHaberl et al. 20061) . Notice, however, that our fits still show 
some variability of the pole temperature correlated to that of 
the orientation angle. In theory, if our model were perfect, it 
would give the appropriate surface temperature distribution to 
explain all variations keeping everything constant except the 
orientation angle. We have given a step forward in that direction 
since our fits are obtained keeping fixed nn, normalization (i.e. 
distance), and the star model (mass, radius, temperature distri- 
bution geometry), but still need to allow for small variations 
of temperature. Another alternative could be that for some un- 
known reason (magnetic field evolution, diffusion) the temper- 
ature and/or size of the hot regions has changed in a timescale 
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of 7 years, but we are not aware of any mechanism able to do 
this in such a short timescale. If this happens to be the case, no 
periodicity would be expected. 

Since we have only explored a few magnetic field config- 
urations an extensive study is very likely to improve the sta- 
tistical quality of the fits presented here. Nevertheless, we are 
confident that some qualitative features are quite robust: i) The 
models that reproduce the observational data have magnetic 
fields confined to the crust and the outer regions, and strong 
toroidal components (about an order of magnitude larger than 
the poloidal one). This naturally leads to large surface tempera- 
ture variations, ii) Models with strong quadrupolar components 
are in general favored. Higher order multipoles (octupole, etc) 
can in principle be added, but a too much patched geometry is 
in conflict with the large pulsation amplitudes, iii) To explain 
the anti-correlation of the hardness ratio we need temperature 
distributions which break the north/south asymmetry. In par- 
ticular, models with a strong quadrupole predict two hot spots 
with different temperatures and sizes and a hot belt that can be 
shifted from the equator 

The emission properties due exclusively to the condensed 
surface model under consideration are the presence of the spec- 
tral edge at about 0.3 keV and the absence of other spectral 
lines expected from gaseous atmospheres with heavy elements. 
We think this part of our results is not as robust as the strong 
anisotropy in the temperature distribution produced by toroidal 
fields and non dipolar components. There are interesting al- 
ternatives that can explain this spectral features equally well, 
but the physical motivation of the condensed surface is that the 
magnetic field is high enough to produce the condensation, if 
the composition is primarily iron or similar heavy elements. 
Strongly magnetized hydrogen atmospheres (or other light ele- 
ments) are still an alternative, but to our knowledge there is no 
fully consistent treatment for arbitrary magnetic field orienta- 
tions and structure. More work in both lines will be certainly 
worth to do in the near future. 

If future observations confirm that this INS (or others) 
is subject to precession with a relatively large wobbling an- 
gle of ^ 20°, combined with a relatively short period of 
~7 yr, it would have very interesting implications. As pointed 
out by Link ll2003h . a large wobbling angle with preces- 
sion timescale of the order of years is in conflict with the 
standard picture of a type II superconducting core. This 
scenario favors either type I superconductivity or a situa- 
tion in which the magnetic field does not penetrate into the 
core and it is mainly confined to the crust and outer re- 
gions. We point out that this was exactly the starting as- 
sumption of our model, a magnetic field living in the crust 
needed to explain the lar^e surface temperatu re anisotropies 
jPer ez-Azorm et al. 2006*: Geppert et al. 2006). Alternatively, 
Jones (2004) suggests that long-term precession leads to the 
conclusion that nuclei and superfluid neutrons do not coexist, 
which also would establish severe constraints on the equation 
of state. At present, the uncertainties in our knowledge of vor- 
tex dynamics in the inner crust can lead to different conclu- 
sions, all of them interesting and based on the observational 
confirmation of long-lived under-damped precession. 



We also point out that using all available constraints from 
the X-ray spectra, the pulsed fraction, the pulse profiles and 
hardness ratios, together with their long-term evolution seems 
a very promising way of constraining the theoretical models. 
We plan to use all the available constraints in the near future 
to relax some of the assumptions made in this work and to ex- 
plore in much more detail the vast parameter space. We are 
also working to study other NSs that we expect to be well de- 
scribed by these models. Preliminary results on RBS1223 are 
promising, since we could find good fits to the very irregular 
light curve and the large pulsed fraction, again with quadrupo- 
lar magnetic fields. Another interesting NS is the pulsar PSR 
J1119-6127, that has a very high pulsed fraction (70 %). This 
is impossible to reproduce with purely dipolar components but 
we found some (still axisymmetric) magnetic field configura- 
tions that can explain this extremely large variability. Results 
about this objects and other isolated neutron stars will be re- 
ported in future work in preparation. In parallel, the issue of 
the magnetic field stability and its diffusion timescale remains 
open. Our final goal is to find a general solution pointing to- 
wards a unified picture in which INSs are just old, cold magne- 
tars whose magnetic fields are a few times smaller than usual, 
maybe because they have decayed during their lifetime. 
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Table 1. XMM-Newton observation information. All observations are performed in Full Frame with the Thin Optical filter 
applied except in Rev. 175 where the Medium filter was used. 



Rev. 


Instrument 


Epoch 


Exposure 


78 


mosl 


2000 Mav 1 ^ 


48 ks 


175 


nn 


2000 Nov. 21-22 


23 ks 


533 




2002 Nov. 6-7 


26 ks 


533 


mosl 


2002 Nov. 6-7 


30 ks 


534 


pn 


2002 Nov. 8-9 


27 ks 


534 


mosl 


2002 Nov 8-9 


32 ks 


622 


mosl 


2003 May 2-3 


24 ks 


711 


mosl 


2003 Oct. 27-28 


14 ks 


815 


pn 


2004 May 22-23 


22 ks 


815 


mosl 


2004 May 22-23 


26 ks 


1086 


pn 


2005 Nov 12-13 


34 ks 



Table 2. The pn and MOS 1 data for all observations are fitted in the 0.18-1.2 keV band with an absorbed BB model including a 
Gaussian absorption line. 



Rev.-Inst. 




kT 






X^/dof 




10-" cm-2 


(eV) 


(eV) 


(eV) 




078-mosl 


0.72!«i 


84 4+'-' 


340+140 

-165 


< 20 


58/53 


175-pn 




83.7+0-^ 


^^^^-105 


< 40 


141/152 


533-pn 


1 20+°" 


85.0!|;i 




32+11 


178/160 


533-mosl 




87.0!|« 


301+^' 


35+25 

JJ_24 


58/50 


534-pn 




85.7+|-| 


308+^5 


351 


159/162 


534-mosl 




86.5+^-3 


320+^0 


25!- 


69/51 


622-mosl 




86.5^2:2 


350+"" 


< 50 


46/48 


711-mosl 


i.80!j:?^ 


903*11 


370!^° 


45+23 


47/47 


815-pn 


1 Q9+0.35 
^•■'^-0.19 


91.4+[-^ 


308+J' 


70+^ 


160/163 


815-mosl 


1 79+0.63 
^- '^-1.05 


91.0^^ 


370!2« 


50!- 


50/52 


1086-pn 


1 oq+0.30 
^•-'^-0.17 


89.6!i;^ 


315+' 


561 


201/170 



Table 3. Joint fits to all EPIC pn and MOS 1 for FF mode observations in the 0. 1 8-1 .2 keV band with absorbed BB plus absorption 
Gaussian line model. The absorption is forced to be the same in all observations. NC indicates that a parameter cannot be 
constrained (i.e. the corresponding model is not required by the data). 



Rev.-Inst. 




KT 




0" line 


-EWune 


^00 / 0300 


^Vd.o.f. 






(eV) 


(eV) 


(eV) 


(eV) 


(km/300 pc) 




175-pn 

533- pn 

534- pn 
815-pn 
1086-pn 




84.45!«:°^ 

S5.451ZI 
85.641- 
91 61+""^ 

^-'•"'--0.04 


NC 

289!|i 
31i:[? 
307+5 
3011 


75/ 


NC 

■xi 0+3.0 

34 4+' ' 
69.8+1^ 
59.3i;^ 


5 5+'-2 
-'•-'-1.0 

5 6+'-5 

J.O_j 4 

5 5+'-' 
5 0+'-^ 
5 1+'-° 

J.l_j2 


902/811 



J.F. Perez-Azorin et al: A self-consistent model of INSs: RX J0720.4-3125. 



13 



Table 4. X-ray spectral analysis. Individual fits to the EPIC-pn and EPIC-MOSl observations in the energy band 0.18-1.2 keV 
for a purely dipolar magnetic field structure (fid = 1.0), ju = 1.34 and O = 12°. We allow to vary the magnetic field intensity at 
the pole (Bp) and the unredshifted temperature at the pole {kTpoie). 



Rev.-Inst. 


nH 








Hp 


/Q.o.r. 




{L\j cm ) 




(,eV) 


(Km/^'-'^' pc J 


( 1 ni2 




78-mosl 


55+°"' 


58.0:j^/ 


118.r3] 


14 1+10.0 


8 4+3 'i 


54/53 


175-pn 


1 06+° " 


52.7+^*2'' 


115 1+*'' 


1 c 9+20.4 
iJ-^_9.4 


10.7!^j 


134/152 


533-pn 


2 15+°'^ 


57 4+^ ' 

' -^-10.8 


111 5+3 3 


in q+11.0 
^^•^-112 


17.4:5-^ 


175/160 


533-mosl 


r» 7Q+0.96 
' ^-0.28 


26 1+3''' 


122.4+2° 


11 3+''-2 


24.9+^°,2 


57/50 


534-pn 


1 Ql+O.ll 

-0.20 




113.3+1* 




21 5+3' 


157/162 


534-niosl 


1 00+°" 

l.UU_Q27 


63.7!«2 


113.9+^'^ 


,74+11.7 
^ '-^-13.4 


14 3+2''-2 


68/51 


622-mosl 


2 60+°''' 


74 8+''-^ 


108.3+3:3 


25.3+?-3 


15 0+2' * 

-10.0 


45/48 


711-mosl 


9 TQ+0.65 
^•-^^-0.31 


54 2+**-' 


123.6+[2 


15.1+^1/ 


25.0+^^ 


50/47 


815-pn 


9 cQ+0.09 
-0.09 


44 1+^-3 
^^•^-1.8 


124.8!^:° 


i3.9:5| 


25.0!°! 


193/163 


815-mosl 


2-4o:°i 


55 1+3-5 
-5.0 


123.8!?:° 


14 9+''-3 
■^^•'-lo.i 


25.0!3-o 


55/52 


1086-pn 


2.68!°:°^ 


49.2:;j 


121.9!°:^ 


16 l^''^ 


25.0!°J 


232/170 



Table 5. Joint fits to EPIC FF mode observations for a purely dipolar field (fid - 1) with /i = 1.34 and O - 12°, forcing the 
hydrogen column density, the normalization constant and the magnetic field intensity to be the same for all observations. 



Rev.-Inst. 


(10^° cm-2) 


S 


(eV) 


R/D30Q 
(km/300 pc) 


Bp 
(10'2 G) 




175-pn 


1 Q^+0.04 
•-^■^-0.02 


20 5+*-' 


120.6^2 


11 3+'-' 
^^•-'-1.8 


14 1+0.6 




533-pn 




29.7![« 


123 3+°-3 








534-pn 




31 3+'-'' 


123.5!°-^ 








815-pn 




56 0+°'' 


132.3!°-5 








1086-pn 




52.2!|J 


131.2+0:3 






1570/819 



Table 6. Individual fits to the EPIC-pn and EPIC-MOS 1 observations in the energy band 0.18-1.2 keV for a quadrupole dominated 
magnetic field (fid = 0.05) with fi = 3.87 and O = 11°. We allow to vary the magnetic field intensity and the temperature at the 
pole. 



Rev.-Inst. 




s 




R/D300 


Bp 


;t^/d.o.f. 




(10^° cm-2) 




(eV) 


(knV300 pc) 


(10'2 G) 




78-mosl 


43+°-3'' 

"•^-^-0.27 


56.0!««-? 


105.6!f; 


12-5!^:; 


j4 2+21.9 


56/53 


175-pn 


89+° '3 

"•°^-0.14 


70 7+''':° 

'"•'-18.0 


104.4!°J 


13.3!^-^ 


39 6+"-* 
-11.7 


138/152 


533-pn 


1 7 1+0.08 
^•'^-0.10 


51 2+3'' 

-6.7 


104 2+°-^ 


16.2+«-? 


<4 4+6.6 


173/160 


533-mosl 


0.83!°:^^ 


Ti 4+36.4 
-"^•^-13.3 


106.4!^3 


15 2+" 


94 Q+20.0 


57/50 


534-pn 


1.65!°;- 


47 4+2.6 


104.1!°:* 


16 4+*^ ° 


55 0+'-° 

-^■^■"-104 


160/162 


534-mosl 


41+°" 

"•^^-0.40 


43 5+''" 

^■^•-'-26.1 


104.9![:* 


14 4+10.0 


43 7+"-3 


69/51 


622-mosl 


Q 44+0.64 
"•^ -0.39 


eg O+108.0 
-'^•■'-49.6 


106.6!|:« 


11.8!^2^' 


94 Q+22.7 
^^■^-19.8 


45/48 


711-mosl 


2 56+°-3' 

^ -0.88 


09 7+6.14 
-'^•'-82 


107.7+[5 


17 2+" ' 

^ '-^-11.9 


55 0+5-° 


50/47 


815-pn 


1 QQ+0.09 
^•^^-0.11 


25.6!!:^ 




15 4+30 


25 0+'-2 


206/163 


815-mosl 


2 00+°-2'' 


36.7!" 




15 4+'-' 
-7.5 


54.0+^0 


54/52 


1086-pn 


1 90+°°3 


25.4!;:2 


109.5!°J 


16 5+3-2 


25.0!°:? 


255/170 
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Table 7. Joint fits to EPIC-pn and EPIC-MOSl FF mode observations for a quadrupole dominated magnetic field (fid = 0.05) 
with = 3.87 and O = 11°, forcing the hydrogen column density, the normalization constant and the magnetic field intensity to 
be the same for all observations. 



Rev.-Inst. 


hh 
(lO^" 


S 


(eV) 


R/D300 
(km/300 pc) 


Bp 

(10'2 G) 




175 
533 
534 
815 
1086 


1 9 Q +0.04 
• -^-0.04 


7Q Q+6.0 

58.6!^:« 
57.6!|« 
26 9+' '' 
31.8!?:' 


104.6!«-^ 
107.5!0-2 

107.5+°':^ 

ii6-3:2j 


12 7+'° 


18-2:1:5 


1192/819 



